ON PERPETUITIES RELATED TO 
THE SIZE-BIASED DISTRIBUTIONS 



Aleksander M. Iksanov 

Abstract. We study perpetuities of a special type related to the size-biased dis- 
tributions. Necessary and sufficient conditions of their existence and uniqueness are 
obtained. A crucial point in proving all results is a close connection between perpetu- 
ities treated in the paper and fixed points of so-called Poisson shot noise transforms. 



Introduction. 

Let V + be the set of all probability measures on the Borel subsets of R + := [0, oo) 
and, for fixed m > 0, := {v G V + : J °° xvidx) — m}. Recall that given 
distribution v G , a distribution v s b is said to be the size-biased distribution 
corresponding to /y, if 

v s b(dx) = rri~ xv(dx). 

Throughout this paper the symbols £(•) and C s b{-) stand for the probability dis- 
tribution and the size-biased distribution of a random variable (rv) in question, 
respectively. 

Consider the distributional equality 

(1) X = AX + B, 

where a random pair (A,B) is independent of an rv X, and "=" means "equality 
of distributions" . In the recent literature it is customary to say that the rv A is a 
perpetuity. Cf. Embrechts, Goldie (1994) for more details. 

In this paper we treat a quite special type of equality (1) in which C{B) = fi, 
C(X) = n s b, and A, B are independent rv's. It is convenient to put B = 77, X = r) s b 
and rewrite (1) as follows 

(2) i] sb = Ar] sb + f). 

Pitman, Yor (2000, p. 35) mention the following problem: "given a distribution of 
A... whether there exists such a distribution of 77". 

In all Propositions stated below we assume that P(^4 = 0) — and keep the 
following notations /1 := C(rj), p := C(A) where 77 and A satisfy (2). First we 
provide necessary and sufficient conditions of the existence of non-zero fj, and reveal 
that when exists this distribution is unique up to the scale. 
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Proposition 1.1 (Existence and uniqueness). 

Assume that ElogA exists, finite or infinite. Then non-zero p exists iff 



(3) 



ElogA < 0. 



Given m > there exists the unique distribution p of mean m. 

In fact the above Proposition is essentially based upon a close relation between 
solutions to (2) and fixed points of so-called Poisson shot-noise transform (see Sec- 
tion 2 for definition and some properties). This is the content of Proposition 2.1. 

Given p = C(A) 6 V + define a random map U p : — ► as follows: 



where 9 = C(X) e P+, C(X s b) = 9 sb , A and X sb are independent rv's and 
"*" stands for the convolution of measures. To be sure that this map is well-defined 
it suffices to verify that U p 9 is of finite mean. To this end consider the measure 
N given by the equality xN(dx) := C(AX s b)(dx). Clearly it is the Levy measure 
of some infinitely divisible (ID) distribution k, say. Since J Q xN(dx) < oo then 
J °° xn(dx) < oo, and now it is obvious that k = V p 9. 

In the proof of Proposition 1.1 it will be shown that /i is a weak limit of iterates 
Up = Up(U p 1_1 ), n E N, and therefore is a fixed point of U p . As it turned out, by 
putting to p some additional moment restrictions one may prove that U p is a strict 
contraction acting on some complete metric space (T, r). Using the Banach Fixed 
Point Theorem allows us to assert that the iterates of V p 9 converge exponentially 
fast to a unique fixed point of U p , for every 9 G T. 

For fixed l<A<2,m>0 consider the set of probability measures 



defined for all 1/1,1/2 £ "^(A): is a metric on 7 ? +(A), and (P^A), ta) is a complete 
metric space. 

In the next Proposition we study the case of p having some finite moments, so 
the function g(x) := EA X is finite (and log-convex) at least on some neighbourhood 
of the origin. 

Proposition 1.2 (Contraction properties and results on moments), a) Let 

(4) g(p) < 1, for some p > , 

and q € (1,2) be any fixed number such that g(q — 1) < 1. Then 

1) Up defined on a complete metric space {V^ l (q),r q ) is a strict contraction, and 

therefore for every 9 £ V+(q) the sequence U™# 7 n = 1, 2, ... converges in r q -metric 



£ sb (U p 9)=U p 9*£(AX sb ) 




In view of Lemma 3.1 by Baringhaus, Griibel (1997) the quantity 
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(hence, weakly) at an exponential rate to a unique fixed point p, ofl] p ; 

2) the rv r\ with C(rj) — p solves (2), Ei] — m and Ei] 1+P < oo. 

b) If there exists non-zero i] satisfying (2) and E?7 P+1 < oo, for some p > 0, then 



It is easily seen that all p's are ID. This follows from the representation of 
non-negative ID distributions due to Steutel (see Sato (1999, Theorem 51.1) for 
a detailed proof). In the next assertion we study the structure of p, as an ID 
distribution more carefully. 

Proposition 1.3 (Infinite divisibility). All non-zero p's are ID with the shift 
and the Levy measure M{dx) — x £(Ar) s b)(dx). Furthermore, p's are compound 
Poisson provided x~ x p(dx) is integrable at the neighbourhood of zero. 

Two other results of the paper deal with the tail behaviour of p. This is changed 
considerably according to whether the rv A may take values greater than 1, or not; 
if not, then whether esssup p = 1 or not. 

Proposition 1.4 (Exponential moments). 

a) pi's have finite exponential moment iff esssup p < 1; 

b) if esssup p < 1 then p's have entire characteristic functions. 

We also give a very simple independent proof of the next result one implication 
of which follows by Proposition 1.4(a). Corollary 4.2 of Goldie, Griibel (1996) 
contains a result concerning general perpetuities in the spirit of the converse part 
of Proposition 1.5. 

Proposition 1.5. Condition "p is concentrated on (0,1]" is necessary and suffi- 
cient to ensure the existence of the solution r\ to (2) which is completely determined 
by its moments. 

Notations and convention. "LT" ("LST" )-Laplacc (Stieltjes) transform, "rv"- 
random variable, "ID"- infinitely divisible, "a. s. "-almost sure(ly), "w.l.o.g." - with- 
out loss of generality; we always take distribution functions to be right-continuous. 



Throughout this Section we assume that all rv's involved live on a common 
probability space (ft, T, P). Let £, £1,^2, •■• be non-negative iidrv's, independent of 
the Poisson flow {t^},? G N with the intensity < A < 00. Given £(£,) fix a Borcl 
measurable function h : (0, 00) — > [0, 00) enjoying the property J °° E[l A h(s)£\ds < 
00. Recall that under the above assumptions random series £i^( r i) converges 

a.s., and its distribution is called a (Poisson) shot noise distribution. The function h 
is said to be the response function. We refer to Vervaat (1979) and Bondesson (1992, 
Section 10) for some additional information regarding the shot noise distributions. 

For a fixed A, consider a Poisson shot noise transform (SNT, in short) Th,\ 
acting on the set (domain) 



EAp < 1. 



Connection to fixed points of shot noise transforms. 




with values in V + as follows 
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Thus the domain Vt is the set of possible distributions for rv £ that would ensure 
the well-definedness of shot noise distribution. By fixed points of the SNT we mean 
non-zero distributional solutions to the equation 

M* = Th,A(M*)' 

where p,* = £(£). We will essentially make use of an equivalent definition of fixed 
points given via the LST ip(s) = e~ sx p*(dx). Namely p* is a fixed point of the 
SNT T h , x iff 



(5) ip(s) — cxp ^A J ((p(sh(u)) — l)du 
In what follows it is assumed that 

(6) the response function h is right-continuous and non-increasing. 

This assumption permits to define the right-continuous and non-increasing gener- 
alized inverse of h given as follows h*~ (z) = inf{w : h(u) < z}, for z < h(0 + ), and 
0, otherwise. Note that fixed points of the SNT for h's concentrated on [0, 1] and 
satisfying (6) have been studied in Iksanov, Jurek (2002). They obtained necessary 
and sufficient conditions of their existence and uniqueness. 

The key ingredient for the proof of all assertions in Section 1 sounds as follows 

Proposition 2.1. a) If n is an rv satisfying (2) with p = C(A) concentrated on 
[a, b] and such that E log A < then p = C{n) is a fixed point of the SNT T^i 
(fi = ii*) with h defined via its generalized inverse 

,6 

(7) h^{x)= J z^ 1 p(dz) 1 x G (a, 6). 

J X 

Therefore h is right- continuous and non-increasing with / °° h(z)dz = 1 and 
/ °° h(z) \ogh(z)dz < 0. 

b) Conversely, if pf is a fixed point of the SNTTh,x with h satisfying (6) and 

/>oo />OC 

(8) A / h(z)dz = I, l h(z) \ogh(z)dz < and lim h(z) = b e (0, oo] 

Jo Jo z ^+° 

then an rv r\ such that C{n) = \i* verifies (2) with an rv A whose distribution 
p{dx) = —\xh^{dx) is concentrated on [a,b]. 

The proofs. 

The first Lemma of this Section is implicit in Athreya (1969, Theorem 1). 

Lemma 3.1. Let <fii(s) and <fi2(s) be the LST's of probability measures of the 
same (finite) mean. If for large enough positive integers n the function ip(s) — 

P 2 ( s ^\ satisfies the inequality ip(s) < Mip(C n s) , where {C n } is a sequence 

s 

of rv's tending to a.s., then <fii(s) = ^(s)- 

The next Lemma is both an existence and uniqueness result concerning fixed 
points of the SNT. In fact, the main part of Proposition 1.1 is merely a combination 
of Lemma 3.2 and Proposition 2.1. 
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Lemma 3.2. Let h satisfies (6), A / °° h(z)dz = 1 and J °° h(z) \ogh{z)dz < 0. 
Then given m G (0, oo) Th,\ has a unique fixed point fi* with J °° x^i*{dx) = m. 

Proof. For fixed to > consider the set of probability measures '■— {v G V£ : 
J °° xv(dx) — to}. Starting with fio = S m , construct the sequence 

Hn ■= Th.AMn-i := T h,\Ho,n = 1,2, ... 

which is trivially well-defined on provided J °° h(z)dz < oo. The corresponding 
LST's tfn\s) = / °° e~ sx [i n (dx), n = 0, 1, ... satisfy equations 

/>OC 

(9) 4 1} ( S ) = e - ms -cxp{-A / (l- ¥ ) n _ 1 («M«)))d«},n=l,2,... 

Jo 

Similarly, for fixed so > 0, let us define yet another sequence {(fin } that satisfies 

(9) for n = 1,2, but 1 — <^ (s) = (1 — e~ s )l{ se [o. So ]}, where 1^ is the indicator 
of set A. 

Let us verify that the weak limit of fi n , as n —> oo, exists and has mean to. As it 
is well-known, this will mean that Th,\ has a fixed point (which is this weak limit) 

In what follows we use some ideas of Durrett, Liggett (1983, the proof of Theorem 
2.7). From (9) one gets, <Pi\s) > <Po\s), i = 1,2 that implies 

^W>AW, n=l,2,.., a > 0,* = 1, 2. 

Thus the monotone and bounded sequence {<p£)},n = 1,2,... has a unique limit 
pW, i = 1,2, with 

(10) limsups-^l - V3 w (s)) < to, i = 1,2. 

For what follows it is essential that i^q 2 ^ (s) > ip^ (s) implies 

(11) - P {2) {s)) < s-\l - ^ (1) («)), s > 0. 

Note ip (1) (s) is the LST of a probability measure ^\ say. Since h(z)dz < oo 
then by dominated convergence it is easily seen that ip^ 1 ' (s) satisfies the fixed point 
equation (5) or equivalently is a (possibly degenerate at 0) fixed point of the 
SNT. It remains to check that G V^ m . 

To this end for n = 0, 1, ... put $„(s) := e s (l - ip^\e~ s )), 

(2) (2) 

*„(a) := -e s (l-^(e- s )+log^(e- s )). In view of assumptions of the Lemma 
n(dz) := —Xzh^(dz) is a probability distribution. Let 6, 6\, 62, ■■■ be independent 
rv's with this distribution. Under these notations one obtains from (9) by change 
of variable and monotonicity of ^ n 



(12) $ n+ i(a) = E$„(s - logfl) - * n («) > E$„( S - log*) - * («), n = 1, 2, . 
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Consider the random walk So = 0, S n = — Y^i=i^°&®i-> n = 1)2,... On iterating 

(12) one gets 

71-1 

$„(«) > E$ (s + S„) -E^$ (s + Si),n= 1, 2, ... 

i=0 

Note that $ (s) := e s (l - exp(- e - ;s ))l{s > logs } and * (s) := 
e s (e" s - (1 - cxp(-e- s )))l{s > logs }. Since Elog(9 t = 

= A J °° /i(z) log h{z)dz < then by the strong law of large numbers S n — ► +oo a.s., 
as n — > oo. Consequently by dominated convergence 

(13) limE$ (s + 5„) = m. 

n—>oc 

Let us verify that ^E'o('S) is directly Riemann integrable (dRi). This together with 
the renewal theorem on the whole line (Feller (1966), Theorem XI. 1, p. 368) will 
imply 

OC 

(14) lim EV$ (s + Si) = 0. 

s^ + oo ^ — ' 

Since e~ s ^o{s) is a decreasing function in s and 

\l/o(s)c?s = J^ s~ 2 (s — (1 — e~ s ))ds < oo (the latter integrand is equivalent to 
a constant near 0) then ^o(s) is indeed dRi. Now according to (13), (14) 

m < liminf lim $„(s) < lirninfs _1 (l - p^(s)) < liminfs _1 (l - (p (1 \s)) 

whence from (10) lim^s -1 ^ — (p^(s)) =mor cquivalently p^ — p* £ "P^ m . 

To prove uniqueness let us assume on the contrary that there exists another LST 
tp(s) with lim s _1 (l — tp(s)) = m that satisfies (5). Choose in Lemma 3.1 

s^+0 

ipi := (p, tf2 :— (f^\ and for each positive integer n, C n := 0\...6 n . Note that 
C n = exp(—S n ) — ► a.s, as n — > oo. By using an elementary inequality 
\er x — e~ y \ < \x — y\,x,y £ R one concludes from (5) ip(s) < Eip(9s), s > 0. 
On iterating the latter inequality n times one gets ip(s) < E^(C„s) which in turn 
implies that Lemma 3.1 does apply. Hence <p = tp^ which completes the proof. □ 

Now we are ready to prove Proposition 2.1. 

Proof of Proposition 2.1. We will prove both parts of the Proposition simultane- 
ously. First of all note that the statement " p is a probability measure" and condition 
(6) with A J Q °° h(z)dz = 1 are equivalent. 

Further suppose that the conditions of part (b) are in force and the SNT Th,x has 
a fixed point p*. Hence <p*(s) — / °° e~ sx p*(dx) satisfies (5), that is, (p*(s) = 

/•OO r-b 

= exp{-A / {l-cp*(sh(u)))du} = exp{\ (1 - ip* {sz))hT {dz)} = 

J0 J a 

— exp{ — A / (1 — ip*{sz))z 1 p(dz)}, 

J a 
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where p(dz) = —Azh^(dz) is a probability measure. The latter follows from (6) 
and (8). In view of Lemma 3.2 condition (8) implies m — xp*(dx) < oo, for 
some m > 0. W.l.o.g. we may and do assume m = 1, and therefore 
lim s-Hl -¥>*(*)) = 1- 

On the other hand, assume that an rv r\ with E77 = 1 satisfies (2) and the 
conditions of part (a) hold. Then the LT ip{s) = Ee~ sr> solves 

/•OO 

tp'(s) = ip(s) / ip (sz)p(dz). 
Jo 

Note that —ip'(s) is the LST of probability measure p a b(dx) :— xp,(dx). By using 
Fubini's Theorem one has log^(s) = J Q S [log (f(u)] du = Jq J Q °° P (uz)p(dz)du = 
f °° z~ x p(dz)(ip(sz) — 1) or equivalently 

/>OC 

<p(s) =cxp{- / (1 - ipisz^z-^idz)}. 
Jo 

Consider the (possibly er-finite) measure v(dz) = (Az) _1 p(dz) and the (right- 
continuous and non-increasing) function V(x) = — J v(dz). Now define the re- 
sponse function h via its generalized inverse h^{x) = V{x). This implies p(dz) — 
—Xzh^(dz) or equivalently (7). 

It remains to verify that ip* (s) = ip(s) which can be achieved in the same manner 
as in the proof of Lemma 3.2 (use Lemma 3.1). The proof is completed. □ 

Proof of Proposition 1.1. Assume that ¥{A = 0} = and there exists non-zero p 
solving (2), but ElogA > 0. By Theorem 1.6(a) of Vervaat (1979) the latter implies 
that C{rj s b) — S c , for some cel. This is possible if c = 0, the case excluded by us. 
A contradiction. 

Assume now that (3) holds and F{A = 0} = 0. Consider the SNT T M with 
the response function h defined via its generalized inverse h*~(z) as in (7). Since 
ElogA < implies L h(z) log h(z)dz < 0, Lemma 3.2 applies that allows us 
to conclude that given m e (0, 00) T^i has a unique fixed point p* with m — 
J °° xp*(dx). It remains to appeal to Proposition 2.1(b). This finishes the proof of 
Proposition 1.1. □ 

To prove Proposition 1.2 we need the following result. 

Lemma 3.3. Let h be a positive measurable function and for some p > 
f °° hP(u)du < 00. If Eg < 00 then E(J2Zi ZMn)) p < 00. 

Proof. Shot noise distributions are ID, and therefore as it is well-known, 

^(Si^i £*M T j)) p < 00 iff Ji x p M(dx) < 00 where M is the Levy measure. Since 

J °° x p M(dx) = Eg / °° h p {u)du < 00, the assertion follows. □ 

Proof of Proposition 1.2(a). 1) Given the rv A with P{A = 0} = and EA^ 1 < 1, 
q G (1,2), define the SNT Th,i with the response function h whose generalized 
inverse h*~ is defined by (7). By Proposition 2.1(a), first, L h(u)du = 1 from 
which one concludes that the restriction of Th,i to the set P+(g) is well-defined 
and for any 9 6 Vmil) Io° x 0^h,i9)(dx) — m; second, J °° h q (u)du < 1 that in 
conjunction with Lemma 3.3 imply x q (Th^8)(dx) < 00. Thus Th,i maps V+(q) 
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into itself. The key observation for what follows is that for any 9 6 V^iq) and 
p = C(A) where A satisfies the above conditions, the map U p is well-defined and 
moreover 

Up = T M 0. 

Formally this can be verified by recording the corresponding LST's. To prove 1), it 
remains to add that according to Lemma 3.5 of Iksanov, Jurek (2002) the restriction 
of Th,i defined on a complete metric space {V^ l {q) 1 r q ) is a strict contraction, and 
apply the Banach Fixed Point Theorem. 

2) The preceding display and Proposition 2.1(b) imply that fi solves (2). Assume 
now (4) holds and at once note that the assertion E77 = m is obvious. We intend 
to verify that (4) implies E?7 1+p < 00. W.l.o.g. set m = 1 and consider two cases: 

1) p £ N: if p = 1 then (2) implies Er] s b = E77/ (1 — EA) < 00 since by the definition 
of the size-biased distribution Eij < 00. Further we proceed by induction. Assume 
it is already known that for some 1 < k < p, fceN 

Er/ 1+fe = E?7 s fe 6 < 00. 

By assumption (4) the function g(x) is log-convex on (0,p) which implies g(l) = 
EA 1 < 1 for all / e (0,p). Therefore E^4 fe+1 < 1. Now an appeal to (2) allows us to 
write formally 

Er, 2+k = E V l+ k = E(A Vsb + V ) 1+k = 

= i( 1+ t ^^E^Et! 1 ^- 1 = 
i=0 ^ ' 

= EA 1+k ETj 1 + k + ^ ( l + ^EA'Erj^Er] 1 ^-* := 

■.= EA 1+k Er] 2+k +t k , 

which yields Erj 2+k = tk/(l — EA 1+k ) < 00 by noting that tk consists of the finite 
number of finite terms. This completes the study of this case. 

2) p £ N: denote by a the fractional part of p then there exists n G N U {0} such 
that p = n + a. Ifn = the assertion follows from the first part of the Proposition 
where in fact was shown that if E^4 Q < 1 then = C(r/) e V+(l + a). Suppose we 
have already verified that for some 0<k<n, keN 

(15) E V 1+k+a = E V k + a < 00. 

Set k a := 1 + k + a. Now (2) implies that one may write formally 

(E^+i))**- 1 = (E^f)*- 1 = 
= (E(Ai lsb +7 1 ) k -) k " < 
< (EA k ") k ^ {E7 1 k " +1 ) k ^ 1 + {Er) k ») k «\ 

the inequality being implied by the triangle inequality in the space L^ a ■ Now the 
latter inequality can be rewritten as follows 



Er? fe Q +i < E^fca/^ _ (EA*-)*- 1 )*- < 00 
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where finiteness is implied by (15) and log-convexity of the function g(x) that 
guarantees 1 — (EA ka ) fc « > 0. A usual inductive argument completes the proof. □ 

Proof of Proposition 1.2(b). While for p > 1 we use an elementary inequality 
(x + y) p > X p + yP, for x, y > 0, to obtain Er] p sb = E(r) + Ar) sb ) p > Er] p + EA P Erf sb > 
EA p Erj P b ; if p e (0,1), then by a variant of Minkowski inequality {Erf^f/P = 
(E(t] + ArisbYf/P > (Er,Pf/P + {EA p ) 1 / p (Erj p b ) 1/p > (EA p ) 1 / p (Erj p b ) 1 / p . Either of 
them implies the desired. The case p = 1 is trivial. □ 

Proof of Proposition 1.3. Clearly, the proof could be done by using the connection 
to the shot noise distributions (Proposition 2.1). However, we choose to use a more 
direct way. Let w be a non-negative ID distribution with the Levy measure N and 
shift 7 > 0. The following representation due to Steutel (1970, p. 86) is well-known 

r-x r-x 

(16) / yu{dy) = / u[0, x - y]yN(dy) + 7 w[0, 4 
Jo Jo+ 

Note that 7 is merely a mass of an atom at zero of the measure R(dx) :— xN(dx). 
On the other hand, each distribution u satisfying (16) with appropriate N and 7 
is ID. 

Putting m :— 1 let us now rewrite equality (2) in terms of distribution functions 
to obtain 

r-x r-x 

(17) H.b[0,x]= yfi(dy)= n[0,x - y]yM(dy), 

Jo Jo 

where the measure M satisfies the relation 

(18) / yM{dy)= [ v sb [0,x/y}p(dy), 

Jo J a 

or equivalently M(dx) = x~ 1 £(Arj s b)(dx). By comparing (18), (17) with (16) one 
immediately concludes that \i is ID with the Levy measure M. Its shift is because 
P{A — 0} = which implies that the measure S(dx) := xM(dx) is atomless at zero. 
To complete the proof recall that finiteness of the Levy measure is necessary and 
sufficient condition for a distribution to be compound Poisson. In our case this 
requirement reduces to the intcgrability of the measure x^ 1 p{dx) at the origin. □ 

Proof of Proposition 1.4- The part (a) follows from Theorem 1.1(b) of Iksanov, 
Jurek (2002) taking into account Proposition 2.1. Turn now to the proof of part 
(b). W.l.o.g. we only consider the case m = 1. By Proposition 1.3 n is ID with the 
Levy measure M(dx) = x' 1 C(Ar] sb )(dx) . By Corollary 25.8 of Sato (1999) 



(19) 



rOO rQO 

Ee srhb = J e sx xfi(dx) < 00 iff J e sx xM{dx) < 00. 



Put c =esssup p € (0, 1). It suffices to verify that if Ee sr,3b < 00 then Ee^/ C - I ' ?sb < 00 
as well. In view of part (a), there exists s = s (c) > such that Eexp(s cr? s fc) < 00. 
Then Eexp(soAri s b) < 00 and consequently e s ° x xM(dx) < 00. It remains to 
use (19) which finishes the proof. □ 
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Proof of Proposition 1.5. Suppose A e (0, 1] a.s. Then KA P < 1, for all p > 0. 
Consequently, by Proposition 1.2 (a) Kr/ P < oo, for all p > 0. In the sequel "— >" 
means "uniquely determines". Let us imagine the chain and prove its validity 
{E?7™} rl(E N — * {^A n } ne j^f — ► C{A) — > £(r/), where is the rv with fixed mean 
to < oo. 

1) {Ery"}„ e N — * {EA™}„ s n via the relation 



2) {EA"}„ eN -> £(A) since A e (0, 1] a.s. (the moments problem). 

3) £(A) — > £(77) which is obvious by Proposition l.l(a)(or 1.2 (a)). 
Conversely, if C(rf) is completely determined by its moments then Er] p < oo, for all 
p > which in view of Theorem 1.2(b) implies KA P < 1, for all p > 0. Clearly, this 
makes A G (0, 1] a.s. and completes the proof. □ 
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